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The possibility of the pion condensation phenomenon in cold and electrically neutral dense bary- 
onic matter is investigated in /3-equilibrium. For simplicity, the consideration is performed in the 
framework of a NJL model with two quark flavors at zero current quark mass and for rather small 
values of the baryon chemical potential, where the diquark condensation might be ignored. Two 
sets of model parameters are used. For the first one, the pion condensed phase with finite baryon 
density is realized. In this phase both electrons and the pion condensate take part in the neutraliza- 
tion of the quark electric charge. For the second set of model parameters, the pion condensation is 
impossible if the neutrality condition is imposed. The behaviour of meson masses vs quark chemical 
potential has been studied in electrically neutral matter. 
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£Nj ' According to a well-known point of view [1|, |2j, |3[ , pionic degrees of freedom and, especially, the pion condensation 
phenomenon might play a significant role in the description of different nuclear matter effects. At present time it is 
widely believed that the dense baryonic matter that might exist inside compact star cores or observed in relativistic 
heavy ion collisions is no more than dense quark matter, obeying an isospin asymmetry. The physics of such a quark 
matter is adequately describes in the framework of QCD with nonzero isospin chemical potential Recently, it 
was shown that a nonzero pion condensate is generated in QCD if /i/ is greater than the pion mass. This result 
I/-} , was obtained in the framework of an effective chiral Lagrangians approach [j] as well as in QCD lattice calculations, 
performed at zero or small values of the baryon chemical potential \ib 0- However, these two approaches are not 
applicable for the description of an isotopically asymmetric matter at moderate baryon density. To overcome the 
problem, it was proposed to study the QCD phase diagram on the basis of Nambu - Jona-Lasinio (NJL)-type models 
0, (see also the reviews 0, @), which contain quarks as microscopic degrees of freedom, in the presence of a 
baryon chemical potential \xb and an isospin /!/ one. In this way the influence of /i^, fij on both the chiral symmetry 
restoration effect lOj and the formation of color superconducting (CSC) dense baryonic matter llj was considered, 
but without taking into account the pion condensation phenomenon. 

Recently, the pionic condensation effect was investigated in some NJL models at nonzero values of [Ib and fij 
k> ! [l2L Il3l Il4| . In particular, it was shown in 01 that quark matter with finite isospin density might exist in two 
different phases. In the first one the baryon density is zero and quarks are gapped, whereas in the second one the 
baryon density is nonzero and quarks are gappless. Note, in [l2,[l3>[lj] the chemical potentials fiB, fJ-i are independent 
external parameters, so the results might be relevant to the physics of the heavy ion collision experiments only, and 
do not describe the real situation inside compact stars. The reason is that matter in the bulk of a compact star should 
be electrically neutral (at least, on average) as well as remain in /3-equilibrium, i.e. all /3-processes that include quarks 
and leptons should go with equal rates in both directions (as a rule, in this case /ii depends o n ufs )- 

In the present paper, we study in the framework of an NJL model, in contrast to |12t I13L ll 1| . the possibility of 
the pion condensation phenomenon in electrically neutral matter with finite baryonic density at zero temperature. 
Moreover, matter in our consideration is required to be in /3-equilibrium. This means that, apart from quark- and 
meson degrees of freedom, it is necessary to take into account charged leptons (electrons only, for simplicity) . Since 
both the pion condensate and electrons have a nonzero electric charge, it is clear that the positive charge of quark 
matter might be compensated in our case by several ways, depending on the competition between electrons and the 
pion condensate. We placed the charge neutrality condition only locally, i.e. suppose that the ground state of matter 
is a uniform phase with zero electric charge density. One should also note that we do not take into consideration 
the CSC effects (see, e.g., EE HE EH), 

so our results are valid for not so large values of the baryon chemical 

potential, say, for [i B < 1200 MeV. 1 

The paper is organized as follows. In Section II the phase structure of the NJL model with zero current quark mass 
is investigated under the requirements of electrical neutrality and /3-equilibrium for the two parameter sets: G = 5.01 



1 The properties of the electrically neutral and /3-equilibrated CSC matter at finite baryon density were investigated in the framework of 
an NJL model in but without taking into account the pion condensation. 



2 



GeV~ 2 , A = 0.65 GeV (set 1) and G = 6.82 GeV" 2 , A = 0.6 GeV (set 2) (G is the model coupling constant, A is the 
three-dimensional cutoff parameter, used in loop integrations). It turns out that for the set I the neutral matter with 
pion condensate is allowed to exist at some values of [i = /is/3, whereas for the set 2 it is forbidden. In Section III 
the mass behaviours of the scalar- and pseudoscalar mesons are considered vs /x for the parameter set 1. 

II. THE MODEL AND ITS PHASE STRUCTURE 

Our investigations are based on the NJL model with two quark flavors. The corresponding Lagrangian has the 
following form 

L q = qi"id v q + G[(qq) 2 + (qtj 5 fq) 2 ] , (1) 

where (i = 1,2,3) are Pauli matrices and, for simplicity, current quark masses are taken zero. Clearly, the 
Lagrangian L q is invariant under transformations of the color SU C (3) and baryon Ub(1) groups as well as under the 
parity transformation P. In addition, this Lagrangian is symmetric with respect to the chiral SU(2)l X SU(2)r group 
(chiral transformations act on the flavor indices of quark fields only) . In particular, it is invariant under the isotopic 
SU(2)i group as well. Moreover, since Q = 1% + B/2 (in the flavor space ^3 = r 3 /2 is the generator of the third isospin 
component, Q = diag(2/3, —1/3) is the generator of the electric charge, and B = diag(l/3, 1/3) is the baryon charge 
generator), the electric charge is conserved too in the NJL model 

Due to the /^-equilibrium requirement, we must incorporate electrons in our consideration. So the full Lagrangian 
of the system looks like 

L = L q + eYid v e, (2) 

where e is the electron spinor field. (We suppose that electrons are free massless particles, for simplicity.) Clearly, 
the Lagrangian @ is well-suited for the description of different processes in the vacuum. To study the properties of 
matter with nonzero baryon- as well as electric charges, we need to modify J5J) as follows 

L = L + f i B N B + f i Q N Q , (3) 

where Nb, Nq are baryon- and electric-charge density operators, correspondingly, and /is, /iq are their chemical 
potentials. 2 Evidently, it holds 

N B = qB^q, N Q = qQj°q - e 7 °e. (4) 

The HQ-term in J3J spoils the vacuum chiral SU(2)l x SU(2)r symmetry of the system. So, at hq ^ the Lagrangian 
L is invariant only under the reduced Ui 3 l(1) X Ui 3 r(1) chiral symmetry group, i.e. the isotopic SU(2)i symmetry 
between u and d quarks is absent in medium. In this case the pion condensation phenomenon might occur. It 
means, without loss of generality, that the ground state expectation value of the form (qi^Tiq) is nonzero, whereas 
(qij 5 T2^q) = (clearly, parity is broken in the ground state of matter with nonzero pion condensate). Another 
characteristic of the ground state of dense matter is the chiral condensate, i.e. the quantity (qq). When it is nonzero, 
the chiral symmetry is spontaneously broken down. In the present paper, in order to establish the phase structure 
of the neutral matter within the framework of the model Q, we restrict ourselves to the consideration of these two 
condensates only. 

The competition between these two condensates is governed by the thermodynamic potential (TDP) which in the 
mean field approximation has the following form (it can be obtained with ease, using, e.g., the technique of (lSjl: 

a x ' 

where the first term in the right hand side is the TDP of free massless electrons. The summation in (J5J runs over all 
quasiparticles (a = u,d,u,d), where 

E u = E A — p,, E a = E A + p,, 

E d = E+- fi, E d = E A +fl, (6) 



2 The Lagrangian can be identically transformed in the following way: L = L + (/ig /3 + fiQ /6)q^ 'q + HQql3l°q — A t Qe7° 
is presented after It is clear from this relation that [iq is just the isospin chemical potential /xr. 
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FIG. 1: The behaviour of the electric chemical potential fig vs quark chemical potential /i = /j.b/3 in the electrically neutral 
matter for set 1 of NJL model parameters. Here fii c ~ 301 MeV, /12c ~ 323.5 MeV. 



and E% = ^(i? ± ) 2 + A 2 , E ± = E±pt Q /2, E = ^Jf- + M 2 , fi = [i B /3 + n Q /6. The factor 3 in front of the summation 
symbol in JSJ indicates the three-fold degeneracy of each quasiparticle in color. Moreover, in order to avoid usual 
ultraviolet divergences, the integration region in J5J) is restricted by a cutoff A, i.e. \p\ < A. First of all, let us fix the 
model parameters as follows: G = 5.01 GeV~ 2 , A = 0.65 GeV (set 1) (Later, another parameter set will be discussed). 
The gap coordinates (M , A ) of the global minimum point of the function f2(M, A) are connected with condensates 
in the following way: 

M = ~2G(qq), A = -2Gi{q 1 5 T 1 q). (7) 

So if Ao is nonzero in the global minimum point (GMP), then the pion condensation phase is realized. Note, the 
quark gap Mo is just the dynamical (constituent) quark mass. From it is possible to obtain the gap equations 

dn(M,A) M f d 3 p ( 9(E+-fl)E+ 9(E A -fi)E- } 

dM ~ 2G J (2nrE\ E+ + E A /' 

n _dn(M,A) _ a__ 6A f d3p ( 6 ( e a-P) , ^Z-m) 



dA 2G J (2tt) 3 I E+ E A 



I Ef E7 V [ ' 



As it was noted in the Introduction, we are going to impose the neutrality constraint locally, i.e. we search for the 
ground state of the system, in which the electric charge density uq = —dQ/d/iQ turns identically into zero. In other 
words, we study the GMP of the function 0(M, A) under the constraint 

= n Q ^^+J -^{0{fi - E+) + 6{il E A ) + 39(E+ - W{E A - (9) 

It is easily seen from the gap equations (jHJ that at fj,Q ^ the global minimum point (GMP) of the TDP |p} might 
take only one of the following three forms in the (M, A)-space: i) (0,0), ii) (M o ,0), and iii) (0, Ao). (In the ground 
state, corresponding to the GMP of the form i), both the chiral- and pion condensates are zero. The solution of the 
type ii) corresponds to the matter phase in which only the chiral condensate is generated. Finally, in the GMP of 
the form iii) the chiral condensate is zero, but the pion condensate is nonzero.) 3 If the neutrality requirement © is 
not taken into account, then the quantities Mq and Ao are functions of the chemical potentials /iq and \xb which are 
independent quantities. In this case the model is relevant for the description of matter with nonzero electric charge 
density. 

However, since we are going to study the electrically (locally) neutral matter, it is necessary to perform the joint 
consideration of the neutrality constraint © and the gap equations JSJl. In this case both the gaps, i.e. the GMP 



3 If hq = 0, then the TDP JSJ depends effectively on the single variable p = VM 2 + A 2 . So, at sufficiently small values of (j,b the global 
minimum of the function Q(M, A) is achieved at all points of some circle in the (M, A)-space. Formally, in this case there is a freedom 
for selecting the GMP. However, since at zero isospin chemical potential, i.e. at fiQ = 0, parity is a conserved quantity in the strongly 
interacting physics, we suppose that in this case the pion condensate is zero, but the chiral one is nonzero (at rather small values of 
/is), i.e. the global minimum is placed by hand in the point of the form ii), and the chirally noninvariant phase is realized (the details 
of the phase structure investigation of the NJL model with zero current quark mass and at (j,q = are presented, e.g., in |2(1 '). 
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FIG. 2: The pion condensate Ao vs /j, — /j.b/3 in the electrically neutral matter for the parameter set 1. Here fii c w 301 MeV, 
ti2c « 323.5 MeV. 

coordinates, and the electric chemical potential /iq are functions of the baryon chemical potential [Ib only (below, the 
analysis is performed in terms of the quark chemical potential /i = /xb/3). Numerical investigations of the equations 
©-© show the following results on the phase structure of the electrically neutral matter for the parameter set 1. 

At sufficiently small quark chemical potential fi < \i\ c w 301 MeV the quantity \iq is equal to zero (see Fig. 1). In 
this case the ground state of the system corresponds to a chirally noninvariant phase. In terms of TDP it means that 
at fj, < Hie the GMP has the form ii) with Mg « 301 MeV. In this phase both the baryon- and the isospin densities 
are zeros, so there is no need to neutralize the quark electric charge. Hence, it is not surprising why /iq = at rather 
small values of /i. 

At /ii c < fi < /i2c ~ 323.5 MeV the GMP has already the form iii), i.e. in this case the pion condensed phase is 
realized. The pion condensate Ao vs \x is depicted in Fig. 2. For this phase of neutral matter the )iq is a nonzero 
negative quantity (see Fig. 1). Nevertheless, the numerical study shows that /xq is a sufficiently small quantity, so 
the relation Ao < p, is fulfilled. As a result, one can see that both u- and d-quasiparticles are gapless in this phase, 
i.e. the quantities E u and Ed from © with A = Ao and M = turn into zero at some energy values, so there is 
no energy cost for creating these quasiparticles. This fact means that at \i\ c < ji < fi2c we have the gapless pion 
condensed phase (GPC) which has a nonzero baryon charge (see also [13 )• I n the GPC phase the electric charge of 
quarks is neutralized both by the electric charge of the electron gas and the charge of the pion condensate. 

Finally, at larger values of the quark chemical potential, i.e. at (i2c < Mi the normal dense baryonic phase is arranged, 
in which both condensates are zero. Hence, in this phase only the electron gas takes part in the neutralization of the 
quark electric charge. Therefore, the absolute values of /iq at fi2c < M are greater, than in the GPC phase (see Fig. 

1) - 

Now, let us consider the phase structure of the electrically neutral matter for another set of model parameters (set 

2) : G = 6.82 GeV~ 2 , A = 0.6 GeV (in this case the dynamical quark mass is approximately 400 MeV in the vacuum, 
i.e. at fiQ, hb — 0.) Numerical analysis shows that for set 2 the phase structure of the model 10 differs qualitatively 
from the set 1 case. Indeed, here at the point /i c as 386.2 MeV we have at once the phase transition from the chirally 
noninvariant phase, which is at fj, < fi c , to the normal dense baryonic phase, at /i c < fi, with zero pion- as well as 
chiral condensates. It turns out that under the neutrality constraint the pion condensation is prohibited in the NJL 
model with set 2 parameters at zero temperature. 4 In contrast, without this constraint the pion condensation is 
allowed to exist in the set 2 NJL model (see [13. IT^ I- 



III. MESON MASSES IN ELECTRICALLY NEUTRAL MATTER 



In the present section the masses of mesons are investigated in the electrically neutral matter. We will follow the 
way, used in plj for studing the particle masses in the color superconducting quark matter. To begin with, let us 
introduce auxiliary bosonic fields 



a(x) = -2G(qq), n a {x) = -2G(gz 7 5 r a? ), 



(10) 



4 This is an important point. Indeed, it was shown in llli Il7l that in electrically neutral matter with rather large baryon density the 
temperature induces a diquark condensation for some range of model parameters, thus tending to the color superconductivity. In a 
similar way, we suppose that at low baryon density the pion condensation can appear in the neutral matter at some temperature interval. 
However, it is a subject of special consideration. 
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where a — 1, 2, 3. In the following we will ignore the influence of electrons on the in-medium meson masses. In terms 
of <j{x) and ^ a { x ) the Lagrangian (with omitted electron part) can be reduced to the form 



L = q 



v • q i - i A*0 • 5 

7 «d„ + ^7 + -^-^37 - cr - ?7 n a T a 



1 

4G 



(7(7 + 7r a 7T„ 



(11) 



(the quantity /i is defined after ©). Starting from (|llf> . it is possible to integrate out the quark fields and obtain the 
effective action of the system in the one-quark loop approximation: 



S c ff(cr, 7r Q ) = - J d A 



where 



D = j"id v 



j d 4 x 






M7° + 


-yr 37 - 



zTr 



sfcx 



ln£> 



«7 ToTa- 



(12) 



(13) 



The Tr-operation in |Q stands for calculating the trace in spinor- (s), flavor- (/), color- (c) as well as four-dimensional 
coordinate- (x) spaces, correspondingly. 

It is clear from JJJ and (|10fl that the coordinates (Mo, Ao) of the global minimum point of the TDP are just the 
ground state expectation values of the a- and 7T]-fields, i.e. Mo = (a(x)), Aq = (iri(x)). 

Let us make the following field shifts in Ijl2|l : — > Mo + cr(x), ni(x) — > Ao + 7ri(x), and thereafter expand the 
effective action up to second order in the meson fields. Differentiating twicely the obtained expression with respect 
to meson fields, it is then possible to get the one-particle irreducible (IPI) Green's functions Txy of the mesons 
(X,Y = a, 7Ti, 7T2, 7Tb) . ( In the present paper we omit these cumbersome calculations, referring to the similar meson 



&) 



mass calculations in |21j.) The results are the following. 

First, let us consider the masses of the a- as well as 7T3-mesons. (Note, tto = tt^.) It turns out that both a- 
and 7To-mesons are not mixed with other particles. Moreover, in the momentum space representation and at zero 
three-momentum, p = 0, we have in the neutral gapless pion condensed phase of matter (/i lc < fi < /i2 C ) : 



r ff(T (po) = r Wo7ro (p ) = -jjj 
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where £ A 



±^2 



(£ A - m) 

2 ,-± _ 



(2tt) 3 



0(4 - fi) 



£aPo + ^Q e+ 



(el-Po) 2 



(s^r (4+Po) 2 



(e A ) 2 



^Q e - £ aPo 



(14) 



+ A 2 ,£ 



( £a - P0 ) 2 ~( £ +) 2 (e A +Po) 2 -(4) 2 

|g| ± ^q/2. Recall, p, = ji + [J.q/6- Eliminating in (|14J) the coupling constant G with 



the help of the A -gap equation (JHJ), we obtain in the GPC phase 

r CTCT (p ) = r W07ro (p ) - (pq - mq)- 



(15) 



Since the zero of a IPI function in the pp-plane defines the mass squared of a particle, it is evident from (I15f) that in 
the GPC phase M l0 = M a — \\iq\ (see Fig. 3, where M CTj7ro are depicted, or Fig. 1 for (iq). 

The expressions for the IPI functions of the a- and ttq mesons in the normal dense quark matter, i.e. at [i2c < 
follow from (|14|l at Ao = 0. The zeros of the Green's functions T a(T (po) and r^^po) in the pp-plane were studied 
numerically in this phase as well. The corresponding mass behaviours vs /i are depicted in Fig. 3 at /i2 C < 

In contrast to the a, 7To-sector, in the sector of 7Ti and TT2 fields the IPI Green's functions (at p = 0) form a nontrivial 
matrix r (po) in the GPC phase, i.e. there is a mixing between m and tt2- 5 Its matrix elements are: 



tZ g Apo) = 6( P 2 - 4A 2 )A( P 2 ), r^(p ) = 6^A( P 2 ), r^(p ) = r^(-po) = 12^0^), 



^GPC 



GPC/ 



^GPCc 



(16) 



where Ao is presented in Fig. 2 and 

A(Po) = 
B{pI) - 



d 3 q 1 




(2tt) 3 1 


,4b 2 -4( £ +)2] 


A 1 


r e +0(4-M) 


(2tt) 3 1 


,4bo-4(4) 2 l 



Abf5-4(£ A ) 2 ] 
£~6>(£ A -jj) 
Abo-4(£ A ) 2 ] 



(17) 



(see also notations after 1)14(11. Note, in the GPC phase Ao < [i and, in addition, \lq < 0. Therefore, the minimal 
value of the quantity £ A in the integrands of (fTTfl is p. As a result, we see that A(p^) and B{j>q) are analytical 



Note, at nonzero current quark mass there is actually a mixing between <r-, 7ri- and TT2 fields in the pion condensed phase. 
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FIG. 3: Scalar- and pseudoscalar meson masses vs /i = /is/3 in the electrically neutral matter for the parameter set 1. Here 
(lie « 301 MeV, M2c « 323.5 MeV. 



functions in the whole pp-plane, except the cut which is at p\ > Afi 2 . Since there is a mixing between m and 7T2 fields, 
the masses of meson modes in this sector are defined by the zeros of the det(r GPC (po)) in the pg-plane, i.e. by the 
equation 

det(r GPC (p )) = 36pg{(pg - 4A 2 1 )A 2 (p 2 l ) - 4B a (pg)} = 0. (18) 

The evident solution of this equation is p\ = 0. It corresponds to a massless meson mode, specified by ttl (see Fig. 3), 
that is actually the Nambu-Goldstone boson. (The appearence of such a mode in the GPC phase, is justified by the 
spontaneous breaking of the initial Ui 3 l(1) X Ut 3 r(1) chiral symmetry down to the abelian subgroup.) The nontrivial 
solution of i|18|) is the zero of the expression in the braces. It corresponds to a massive meson mode, denoted by tth 
(see Fig. 3). (Clearly, its mass, M 1TH1 lies in the interval 2Ao < M„ H < 2ft.) 

In the dense normal quark matter (NQM) phase, i.e. at at fj,2 C < fi, it is convinient to use the charged fields 
7T±0) = (iri(x) ± i7r 2 (x))/V2. Then in the NQM phase the matrix r NO - M (p ) of 1PI Green's functions of the tt± 
mesons looks like: 



r»T- foO - r^ M + (-po) = ^ - J, { a 2 - ^ - 2f» - („ + !£ ) 



2 



(MQ + Po 



■In 



4A 2 -(/XQ+p ) 2 
(2 M + MQ /3) 2 -p 2 



}, r^ M bo) = r^ M bo) = o. (19) 



The numerical investigation of the zeros of the quantity det(F (po)) shows the presence of two pionic massive 
modes, which might be identified in NQM phase with 7r± -mesons (see Fig. 3). Evidently, in this phase the mass 
splitting between 7r±-mesons is due to the isospin asymmetry that is generated by nonzero /iq. 

Finally, in Fig. 3 the masses of mesons in the chirally noninvariant phase, i.e. at fi < fj,\ c , are presented as well. In 
this phase fig = 0, and only the chiral gap Mq ^ 0. It is well-known that in this case M a = 2Mq w 602 MeV, and 
the three 7r-mesons are massless Nambu-Goldstone bosons. 



IV. SUMMARY 



In the present paper we have studied the properties of electrically neutral and /3-equlibrated cold matter with finite 
baryonic density. The problem is inspired by the physics of compact stars. For simplicity, the consideration was done 
in the framework of a NJL model with zero current quark mass. We have found that for the set 1 of model parameters 
(see the end of Introduction) there are three different phases, including the one with a pion condensate, of neutral 
matter. In contrast, for the parameter set 2 the pion condensation in the neutral matter is forbidden. Moreover, we 
have studied the behavior of meson masses vs quark chemical potential in the case of the parameter set 1 (see Fig. 
3). Since the electric neutrality of the system is realized together with an isospin asymmetry between quarks, it turns 
out that the masses of 7r-mesons are splitted at nonzero baryon density. 
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